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5h 1 Abstract 

a 

■ The method of lowest-order constrained variational, which predicts reason- 

ably the nuclear matter semi-empirical data is used to calculate the equation 
of state of beta-stable matter at finite temperature. The Reid soft-core with 

(N 

. and without the N-A interactions which fits the N-N scattering data as well 

m , 

|Q . as the UV14 potential plus the three-nucleon interaction are considered in 

the nuclear many-body Hamiltonian. The electron and muon are treated 
relativistically in the total Hamiltonian at given temperature, to make the 

. 

fluid electrically neutral and stable against beta decay. The calculation is 

Ph 

performed for a wide range of baryon density and temperature which are of 
interest in the astrophysics. The free energy, entropy, proton abundance, etc. 
of nuclear beta-stable matter are calculated. It is shown that by increas- 
ing the temperature, the maximum proton abundance is pushed to the lower 
density while the maximum itself increases as we increase the temperature. 
The proton fraction is not enough to see any gas-liquid phase transition. Fi- 
nally we get an overall agreement with other many-body techniques, which 

are available only at zero temperature. 
PACS No. 21.65.+f, 26.60,+c, 64.70.-p 
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I. INTRODUCTION 



The prediction of equation of state of hot and dense nuclear matter has been a subject 
of growing interest in nuclear physics and astrophysics. Such matter under such unusual 
conditions, beside of various nuclear physics applications, it is likely to be produced through 
stellar collapse, supernova explosion, neutron stars etc. [1]. 

On the other hand, it is believed that at extremely high density, deconfmement can take 
place and the hadronic matter may make a transition to the quark matter. So the properties 
of transition points depend on the equation of state of hadronic matter. 

It is well known that the presence of leptons is also crucial since the star matter have 
to be equilibrated against the weak leptonic-decay. So the increase of charged leptons with 
certain baryonic density will finally turns to a large fraction of neutron density. However, 
other negatively charged baryons such as S~ hypron may also become energetically more 
favorable [2]. 

In order to study such systems one needs a good many-body techniques to provide an 
accurate description of the equation of state, entropy and other thermodynamics variables 
of hot nucleonic matter. 

In a series of papers the lowest order constrained variational (LOCV) method has been 
developed [3-8] for calculating the properties of homogeneous nuclear fluids with realistic 
nucleon-nucleon interaction. This approach has been further generalized to include more 
sophisticated interactions such as the V u [8], the AVu and the new argonne AV\% [8] as 
well as the Reid [3] and A-Reid [3] potentials. For a wide range of models our LOCV 
calculation agrees well with the results of fermion hypernetted chain (FHNC) calculations 
where these have performed and for a number of central potentials there is agreement with 
the essentially exact numerical solutions obtained by Monte Carlo technique [3] . Despite this 
agreement for model problems, there has been some dispute about convergence of LOCV 
results in calculation employing realistic nucleon-nucleon interactions which are strongly 
spin-dependent and which, in particular, contain a sizeable tensor force. This argument 
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was tested by us by calculating the energy of the three-body cluster contribution in nuclear 
matter and the normalization integral < ip \ ip > both at zero and finite temperature [4,7]. 
It was shown that < ip \ ip > is normalized within one percent and the three-body cluster 
energy is less than one MeV for kf < 1.6/m -1 . Our LOCV calculation is a fully self- 
consistent technique and is capable of using the well-defined phenomenological potentials 
such as A-Reid (the modified Reid potential with an allowance of A(1234) degree of freedom, 
see Modarres and Irvine [3]) potential. The A state is most important configuration which 
modifies the nuclear force and it might be at the origin of understanding of three-body forces 
[9]. The results suggest that the LOCV method reasonably describes the nucleonic- matter 
properties at zero and finite temperature. 

On the other hand, our recent calculation at zero temperature [8] with V 14 and V i8 poten- 
tials show the same behavior and a very good agreement was found with more sophisticated 
calculation such as the fermion hypernetted chain method (FHNC) [10]. 

With respect to the above arguments, in this work we shall attempt to calculate the 
properties of beta-stable matter at finite temperature by using the Reid and A-Reid potential 
and investigate the behavior of proton abundance with the temperature and density of 
baryonic matter. 

Brueckner theory [11] and variational FHNC [12] have been also applied to beta-stable 
matter but only for zero temperature. So we can also compare our result with above calcu- 
lations. 

The paper is planed in the following steps: The beta-stable matter free energy is ex- 
plained in section II. Section III is devoted to a short description of the lowest order con- 
strained variational method. Finally in section IV we present the results and discussion. 

II. THE BETA-STABLE MATTER FREE-ENERGY 

The free energy of beta-stable matter is written as the sum of free energy of the leptonic 
(F L ) and the baryonic (F B ) parts : 
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F = F L + F B (1) 

The nucleons are assumed to interact through one of the realistic N-N interactions i.e. Reid, 
A-Reid and Vi 4 potentials. The requirement of charge neutrality implies that we can ignore 
the electromagnetic interactions and the week interactions are neglected. 

The total baryon number density n B is the sum of the proton and neutron number 
densities 

n B = n p + n n (2) 
and the condition of electrical neutrality requires: 

n p = n e + (3) 

(here we ignore r leptons because of its large rest mass respect to two other leptons.) 

The leptons form two highly relativistic fermi seas. The contribution to the energy per 
baryon from these fermi seas is, 

E L = (Qnsr 1 ]T YlHk) + m iC 2 \fi{k) (4) 

i=e,fi k,cr 



where 



£.(£;) = [h 2 k 2 c 2 + m 2 c 4 f 2 - m iC 2 (5) 



and 



f l (k) = [exp((e t (k)- l 2 t )P) + l}- 1 (6) 

is the familiar Fermi-Dirac distribution with f3 = (K is the Boltzman factor). ^ are the 
chemical potentials of the ith pedes of particle. 

For zero temperature, E L can be calculated analytically and it take the following form: 

E L (T = 0)= £ -^J^[a; l (l + a; l 2 ) 1 / 2 (2x l 2 + l)-sinh- 1 x l ] (7) 

where 



x, = ^ (8) 
rriiC 

and k e and /c M are the electron and muon Fermi momenta respectively. The chemical po- 
tentials ^ (or Fermi seas in case of T = 0) are related through condition of beta stability 
i.e. 

fin ~ Hp = He = /V (9) 

For T = the second equality in above equation implies that 

m e c 2 (l + ^) 1 / 2 =77 V c 2 (l + 4) 1 /2 (10) 

while for T ^ the chemical potentials are fixed by the various particle number densities 
i.e. 
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(n _1 )E/« (") 

and the second equality in equation (9) should be fulfilled numerically. It is clear that as 
far as fi e is less than the rest-mass of muon we will not have any muon in the matter. 

The first equality in equation (9) is satisfied by minimizing the total free energy respect 
to the proton number density n p . 

Then the leptonic free energy per baryon can be written as 

F L = E L - TS L (12) 

where 

Sl = S e + S„ (13) 

and 

Si = Kinns)- 1 £(1 - fi(k)) ln(l - + f t {k) HMQ) (14) 

fc,cr 

Finally we write the baryonic free energy as 



Fb — Eb — TSb 



(15) 



where 

S B = S p + S n (16) 

with definitions of Si from equation (14) and Fermi-Dirac distribution of equation (6), but 
h 2 k 2 

with Si = 2—k- m* are the proton and neutron effective masses and they will be treated 

i 

variationally. The baryonic internal energy is 

E B = T B + E B 4B (17) 

where the kinetic energy part has the following form (the baryons are treated non- 
relativistically), 

i=n,p 13 1 

= ^n B )- 1 Y,Y,[^ + m Am if T^O (18) 

k,a i=p,n * 

and k n and k p are the familiar neutron and proton fermi momenta, respectively. The many- 
body energy term E B 4B will be discussed in the next section. 



III. THE LOCV FORMALISM AND E^ B CALCULATION 

The calculation follows exactly that of asymmetric nuclear matter calculation of reference 
[6]. We consider an ideal Fermi gas type wave function for the single particle states. Then 
using variational techniques, we can write the wave function of interacting system as, 

$ = T r $ r (19) 

where 

FT = SHf(ij) (20) 

i>j 
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The Jastrow correlation functions f(ij) are operators that act on spin, isospin and relative 
coordinate variables of particles % and j and S is a symmetrizing operator which is neces- 
sary, since the f(ij) do not commute. The unitariness of jF r 's usually cause the problem 
of nonorthogonality of different states. But since at low temperature only the one quasi 
particle-type states are important and these states have different total momentum, they can 
be considered to be orthogonal. 

Now, using the above trial wave function, we can construct a cluster expansion for the 
expectation value of the following Hamiltonian, 

^ = E-|^ + £^ (21) 

i zm i<j 

where Vij is a two-nucleon potential that fit the nucleon-nucleon scattering data and deuteron 
properties. In this work we will mainly focus on the Reid [3] and A-Reid [3] interactions. 
Since, regarding our previous works, the A-Reid interaction can reasonably reproduce the 
nuclear matter properties [3,8]. We will also use the UV U plus density dependence three- 
body potential in order to compare our result with others [12] at zero temperature. 

In the cluster expansion series, we keep only the first two terms of energy functional: 



E B ([f}) = 4 =E l + E 2 (22) 



1 < V 


H 


ip > 


A <ip 


ip > 



Ei is independent of f(ij) and is simply the Fermi-gas kinetic energy of baryons T B which 
already has been discussed in equation (18). While the two-body energy E 2 is written as 

E 2 = ^ E < V I V ( 12 ) I V ( 23 ) 

where 

V(12) = -^[/(12) , [V? 2 , /(12)]] + /(12M12)/(12) (24) 

and the two-body antisymmetrized matrix element < ij \ V | ij > a is take with respect to 
the single-particle wave function $ r . 

We minimize the two-body energy E 2 with respect to the variations in functional f(ij) 
but subject to normalization constraint [3,8], 



where 



with (i = p, n) 



i^< U '|^ 2 (12)-/ 2 (12)|u> a =0 (25) 



A 



/i = [1 ( — — ) 2 1~2 n — n and p — p channel 

y 2 y n B 

= 1 n — p channel (26) 



7,(r) = ^p//,(A:)Jo(A:r)rfk (27) 

The above constraint forces the cluster series to converge very rapidly and we can approxi- 
mate E 2 by E^ B i.e. the many-body contribution to the internal energy. 

The detail about the Euler-Lagrange differential equations, the boundary conditions and 
the channel break down of E 2 have been discussed in references [3,6]. 



IV. RESULTS AND DISCUSSION 

The equilibrium configuration of beta-stable matter is obtained at each total baryon 
number density % and temperature T by minimizing the total free energy F with respect 
to the two-body correlation functions, effective masses and proton abundance n p /nB subject 
to the constraints of equations (2), (3), (9) and (25). 

In figure 1 we plot our calculated free energy per baryon number at different tempera- 
tures for the Reid and A-Reid potentials. This figure can be compared with our previous 
calculation for pure neutron matter [5] (figure 2). It is seen that beta-stable matter has 
softer equation of state than pure neutron matter and as we expected the free energy per 
baryon decreases with increasing temperature at constant density. But with increasing the 
density the free energy shows an increasing trend for given temperature. Inclusion of isobar, 
increases the free energy and makes it more density dependent. 

The relative proton abundance, y p — for the Reid and the A-Reid potentials at 
different temperatures is plotted in figure 3. It is seen that yp increases at low baryon 

8 



density as we increase the temperature and the maximum value of yp is pushed to the lower 
baryon density by an increase in the temperature. This shows that at high temperatures the 
beta-stable matter with higher proton abundance tends to lower its baryon density. This 
effect is much larger in case of the A-Reid interaction. 

The proton abundance peak is due to this fact that the system takes the advantage of 
the spin-triplet tensor components in the Reid (or A-Reid) interaction and by maximizing 
the proton density at given baryon density makes the free energy lower. Since the A-Reid 
interaction has more repulsive components, like l So — > 1 D 5 transition, the peak is happened 
in lower baryon density. 

In figure 4 the proton, electron and muon number density are plotted against baryon 
number density at KT = 5 and 20 MeV. We find that n e ,n^ and n p increase by increasing 
temperature at low baryonic density while their maximum values become lower as we increase 
the temperature. The entropy of proton, electron and muon are plotted versus baryonic 
density in figures 5 and 6 for the Reid and A-Reid potentials at different temperatures. 
It is seen that the entropy of proton and electron increase very rapidly with increasing 
the temperature at low densities while for the muon's entropy we find a peak for given 
temperature. The baryon density in which the muon's entropy becomes maximum goes to 
lower densities as we increase the temperature. 

In figure 7 we compare our beta-stable results for y p with those of Wiringa et al. (WFF) 
[12] and Baldo et al. (BHF) [11] for UVu and AV U potentials at zero temperature. We 
get quite good agreement with them especially at low baryonic densities. For % > .5 the 
BHF calculation increases while our results and WFF decreases. Our results with Reid and 
A-Reid potentials at zero temperature are also given for comparison. 

We would conclude by pointing out that we have carried out a calculation of beta-stable 
matter using techniques, which we believe to be more reliable than Brueckner calculation 
especially at higher densities {jib > 0.2/m~ 3 ) and it is comparable with the variational 
hypernetted chain method. It was seen that the extra degree of freedom offerd by beta 
decay leads to a slight softening of the equation of state, as one would expect. We found 
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that, the maximum values of proton abundance which are predicted by the Reid and A-Reid 
potentials are much smaller than UV U and AV U potentials and they happened in much 
lower baryonic densities. Increasing the temperature of the beta-stable matter increases 
the proton abundance in low densities. This effect becomes more stronger when A-Reid 
interaction is considered. In general, the proton fraction is not enough to se the gas-liquid 
phase transition. It also dose not reach the critical value of about 15 percent needed for the 
occurrence of direct Urea Processes which believed to be responsible for a fast neutron star 
cooling [1]. However it was shown that the temperature may play a great role to reach to 
higher percentage of proton abundance and lower baryonic density. 

Finally the equation of state of beta-stable is more soften than neutron matter and 
the effect of temperature may influence the behavior of dense matter and the structure of 
neutron star. 
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FIGURES 

FIG. 1. Free energy of beta-stable matter (without baryon rest masses) (MeV) versus baryon 
density (fm~ 3 ) at various temperatures (MeV). Dashed curve and full curves are for the A-Reid 
and Reid potentials, respectively. 

FIG. 2. Same as figure 1 but for neutron matter. 

FIG. 3. The proton abundance for the A-Reid (dashed curve) and Reid (full curve ) potentials 
at various temperatures (MeV). 

FIG. 4. Particle densities (fm~ 3 ) of e, fi and p for the A-Reid (dashed curve) and Reid (full 
curve) potentials at KT = 5 and 20 MeV. 

FIG. 5. The entropy of proton (Mev -1 ) in the beta-stable matter for the A-Reid (dashed curve) 
and Reid (full curve ) potentials at various temperatures (MeV). 

FIG. 6. As figure 5 but for e, fi. 

FIG. 7. The comparison of proton abundance y p versus baryon density (/m~ 3 ) at zero tem- 
perature with those of references [11,12] for te UV±4, AV14, Reid and A-Reid potentials. 
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